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Abstract 



3 ■ 

^ ■ We propose that supernova relic neutrino (SRN) observation can be used to set constraints on 

r**"- . the neutrino decay models. Because of the long distance scale from cosmological supernovae to 

the Earth, SRN have possibility to provide much stronger limit than the present one obtained 



from solar neutrino observation. Since the currently available data are only the upper limit on 
the flux integrated over E De > 19.3 MeV, the decay models on which we can set constraints is 
\ quite restricted; they must satisfy specific conditions such that the daughter neutrinos are active 



species, the neutrino mass spectrum is quasi-degenerate, and the neutrino mass hierarchy is normal. 
Our numerical calculation clearly indicates that the neutrino decay model with (Tv/m^Ts/m) ^$ 



a.: 

(10 10 , 10 10 ) [s/eV], where T{ represents the lifetime of mass eigenstates Vi, appears to give the SRN 
flux that is larger than the current upper limit. However, since the theoretical SRN prediction 

•rH . 

contains many uncertainties concerning a supernova rate in the universe or simulation of supernova 
explosions, we cannot conclude that there exists the excluded parameter region of the neutrino 
lifetime. In the near future, further reduced upper limit is actually expected, and it will provide 
more severe constraints on the neutrino decay models. 
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I. INTRODUCTION 



ic Q, 



A number of ground-based experiments, which observed atmospheric |l|, solar 
and reactor neutrinos 5], have revealed the nonzero neutrino masses and flavor mixings, 
i.e., properties beyond the standard model of the particle physics. Fortunately, our current 
knowledge of the neutrino mass differences and mixing angles further enables us to consider 
:'ar more exotic properties of the neutrino, such as a nonzero magnetic moment (see Refs. 

in 

a, |j| and references therein) and neutrino decay. In this Letter we show that neutrino decay 
of a particular type may be ruled out or severely constrained by the current and the future 
supernova relic neutrino observation at the Super-Kamiokande (SK) detector j^]. 

We consider two-body neutrino decays such as z/j — > Vj + X, where z/j are neutrino mass 
eigenstates and X denotes a very light or massless particle, e.g., a Majoron. The strongest 
limits on this decay modes are obtained from the solar neutrino observation by SK J3] (see 
also Ref. ^ was argued that the limit is obtained primarily by the nondistortion of the 

solar neutrino spectrum, and the potentially competing distortions caused by oscillations as 
well as the appearance of active daughter neutrinos were also taken into account. However, 
owing to the restricted distance scale to the Sun, this limit is very weak, r/m > 10~ 4 s/eV, 
and therefore, the possibility of the other astrophysical neutrino decay via the same modes 
cannot be eliminated. In fact, detectability of the decay of neutrinos from the high-energy 



astrophysical sources was discussed 



and it has been concluded that it should be visible 
by future km-scale detectors such as IceCube, since the neutrino decay strongly alter the 
flavor ratios from the standard one, (p Ue : (f) u : <J) Ut = 1:1:1, expected from oscillations 
alone. 

Our strategy is basically the same as that of the previous studies jj], [ll|, i.e., the en- 
hancement of the electron neutrino events due to the decay is investigated. As a source of 
neutrinos, we consider supernovae. Two observational results concerning supernova neutri- 



nos exist; one is the well-known neutrino burst from SN 1987A 
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3, 



and the other is the 



recent upper limit on the flux of supernova relic neutrinos (SRN), which is the accumulation 



of neutrinos from all the past supernovae, by SK 



} Original discussions concerning the 



SN 1987A signal have been already given, including the effect of the neutrino decay as well 



1 From this point on, we consider only v e at detection because this kind of flavor is most efficiently detected 
at SK via v e p — > e + n reaction. 
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as the pure flavor mixing, in the literatures 



Q (see also Ref. Q 



1J, |15[ (see also Ref. |l6j|). In this Letter, we 



use the latter one (SRN) for obtaining implications for the neutrino decay. The advantage 
of this approach compared to the former ones is that the neutrinos must transit over very 
long distance scale from the cosmological supernovae to the Earth, and much longer life- 
times would be probed in principle. The current SRN upper limit is only a factor three 
larger than theoretical predictions by Ando et al. 17, Q] (hereafter AST), which adopted 



neutrino oscillations using experimentally inferred parameters. Therefore, if some neutrino 
decay model predicts the SRN flux which is three times larger than the AST prediction, 
then it is excluded. 

This Letter is organized as follows. In Section |H1 we present formulation for the SRN cal- 
culation and adopted models. In particular the detailed discussion concerning the supernova 
model and the neutrino decay model are given in Section III Al and III Bl respectively. How 
the oscillation and decay changes the neutrino spectrum and flux is qualitatively illustrated 
in Section Mil and the numerically calculated results are given in Section IIVI In Section 
fV] we discuss the various possibilities of the neutrino decay, uncertainties of the adopted 
models, and future prospects. 

II. FORMULATION AND MODELS 

The SRN u e flux is calculated by 

cj Q R^(z)U,XE'^z){l + z)—dz } (1) 



dE p 

where E' v — (1 + z)E u , Rs^s(z) is a supernova rate per comoving volume at redshift z, 
and z max is the redshift when the gravitational collapses began (we assume it to be 5). 
J\f u (E' u , z) = dN u (E' u ) I ' dE' v is the number of emitted neutrinos per unit energy range by one 
supernova at redshift z. As the supernova rate, we use the most reasonable model to date, 
which is based on the rest-frame UV observation of star formation history in the universe 
by the Hubble Space Telescope jl^ . and the model was also used in AST as "SN1." In 
this model, the supernova rate exponentially increases with z, peaks around z ~ 1.5, and 
exponentially decreases in further high-z region. 
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A. Supernova model 



As original neutrino spectra N„{E,z), which is not the same as M v owing to the neu- 
trino oscillation and decay, we adopt the result of a numerical simulation by the Lawrence 
Livermore group j^j], which is the only group that is successful for calculating neutrino 
luminosities during the entire burst (~ 10 sec). The average energies are different between 
flavors, such as (E Ue ) ~ 11 MeV, (E De ) ~ 16 MeV, (E Vx ) ~ 22 MeV, where u x represents 
non-electron neutrinos and antineutrinos. This hierarchy of the average energies is explained 
as follows. Since v x interact with matter only through the neutral-current reactions in su- 
pernova, they are weakly coupled with matter compared to v e and v e . Thus the neutrino 
sphere of v x is deeper in the core than that of u e and u e , which leads to higher temperatures 
for u x . The difference between u e and v e comes from the fact that the core is neutron-rich 
and v e couple with matter more strongly, through v e n <-> e~p reaction. 

Our calculations presented in Section |IV| are strongly sensitive to the adopted supernova 
model, or in particular, the average energy difference between u e and v x . Recently, the 
Livermore calculation is criticized since it lacks the relevant neutrino processes such as 
neutrino bremsstrahlung and neutrino-nucleon scattering with nucleon recoils, which are 
considered to make the mean energy difference between flavors less prominent; it has been 
actually confirmed by the recent simulations (e.g., Ref. [sij]). However, we cannot adopt 
recent models, even though they include all the relevant neutrino microphysics. This is 
because the SRN calculation definitely requires the time-integrated neutrino spectrum during 
the entire burst, whereas all the recent supernova simulations terminate at ~ 0.5 sec after 
the core bounce. Since the Livermore group alone is successful to simulate the supernova 
explosion and calculate the neutrino luminosity during the entire burst, we use their result 
as a reference model. Again, we should note that our discussions from this point on heavily 
relies on the adopted supernova model. 

B. Neutrino decay model 

In this Letter, we consider only the so-called "invisible" decays, i.e., decays into possibly 
detectable neutrinos plus truly invisible particles, e.g., light scalar or pseudoscalar bosons. 
The best limit on the lifetime with this mode is obtained from solar neutrino observations 
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and is r/m ^ 10~ 4 s/eV |9|, which is too small to set relevant constraints on discussions 
below. We do not consider other modes such as radiative two-body decay since they are 
experimentally constrained to have very long life times (e.g., Ref. |22j|). 

We note that our approach is powerful only when the decay model satisfies specific con- 
ditions such that (i) the daughter neutrinos are active species, (ii) the neutrino mass spec- 
trum is quasi-degenerate (mi m ni2 ~ m.3), and (iii) the neutrino mass hierarchy is normal 
(mi < ma), not inverted (mi > m^). This is because at present, only the upper limit of 
the SRN flux is obtained, and therefore, the decay model which does not give large flux at 
detection energy range (Ep e > 19.3 MeV) cannot be satisfactorily constrained. All these 
conditions (i)-(iii) must be satisfied to obtain severe constraints on the neutrino lifetime, 
because (i) if the daughter neutrinos are sterile species, the SRN flux decreases compared 
with the model without the neutrino decay; (ii) if the neutrino mass spectrum is strongly 
hierarchical, then the daughter neutrino energy is degraded compared to its parent and the 
predicted SRN flux at high energy range would not be large; and (iii) in the case of inverted 
hierarchy, v e most tightly couple to the heaviest mass eigenstates, which decay into lighter 
states, and it also reduces the SRN flux. 

For a while, we assume that the conditions (i)-(iii) are satisfied; all the other possibil- 
ities are addressed in detail in Section IV Al As discussed in Section IV CI we believe that 
future observational development would enable far more general and model-independent 
discussions, which are not restricted by the above conditions. 

III. NEUTRINO OSCILLATION AND DECAY 

Effects of neutrino oscillation and decay are included in Afp e (E, z). We address the prob- 
lem of flavor mixing by the pure supernova matter effect, first for antineutrinos and second 
for neutrinos, and then we discuss the neutrino decay. 

The state of v e produced at deep in the core is coincident with the lightest mass eigenstate 
owing to large matter potential. This state propagates to the supernova surface without 
being influenced by level crossings between different mass eigenstates (it is said that there are 
no resonance points). Thus, v e at production becomes i>\ at the stellar surface (Ap 1 = Af®J 
and the number of v e there is given by 

Mp e (E, z) ~ cos 2 BqN Di (E, z) + sin 2 9 Q M i2 (E, z) 
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= cos 2 6 & Af® e (E, z) + sin 2 6 & ^(E, z), 

(2) 

where 9g, is the mixing angle inferred from the solar neutrino observations (cos 2 # ~ 0.7) 
|2,0,U, [4, A/"° represents the neutrino spectrum at production. In the above expression (J2J), 
we used the fact, sin 2 26 l atm = 1, and assumed 6*13 = 0. These are justified by the atmospheric 
and reactor neutrino experiments (2^ . 

The situations changes dramatically for neutrino sector. As the case for antineutrinos, v e 
are produced as mass eigenstates owing to large matter potential, however, the difference is 
that the produced v e coincide with the heaviest state z/ 3 . Since in vacuum u e most strongly 
couples to the lightest state v\, there must be two level crossings (or resonance points) 
between different mass eigenstates during the propagation through supernova envelope; each 
is labeled by H- and L-resonance, corresponding to whether the density of the resonance point 



is higher or lower (see, e.g., Ref. 24j for details). It is well-known that at L-resonance the 
mass eigenstate does not flip (adiabatic resonance) for LMA solution to the solar neutrino 
problem. However, the adiabaticity of the H-resonance becomes larger than unity when 
the parameter #13 is sufficiently large. Instead, we parameterize the flip probability at the 
H-resonance by Ph, i.e., if the resonance is adiabatic (nonadiabatic), Ph = 0(1). Thus, at 
the stellar surface, the neutrino spectra of mass eigenstates is given by 

M vl {E,z) = K x (E,z), (3) 
K 2 (E,z) = P H N° u SE ) z) + {l-P H )Nl{E ) z) ) (4) 
M^{E,z) = {l-P H )N Q Ve {E,z) + P H Nl{E,z). (5) 

If we include the neutrino decay, the expected v e flux from each supernova changes 
drastically. Before giving a detailed discussion, we first place some simplifying assumptions. 
Instead of the lifetime, we define "decay redshift" zf of the mass eigenstate fiipi); if the 
source redshift z is larger than the decay redshift zf, all the neutrinos Vi[pi) decay, on the 
other hand if z < zf, v^ipi) completely survive. We consider the decaying mode z/ 3 (^ 3 ) — ► i>\ 
and ^2(^2) - > v\i and z% and z$ are taken to be two free parameters. The other case that 
one of them is stable can be realized if we take z d > z max . With these assumptions and 
parameterization, the neutrino spectrum which is emitted by the source at redshift z can 
be obtained. First, we consider the decay mode Pi 9j +X, in which the neutrino helicity 



conserves. The v e spectrum is given by 



Nu K {E, z) = cos 2 6 & Af Pl (E,z) + sin 2 6 Q M 92 (E, z) 



COS dr. 







'Af° 1 (E,z)+N° 2 (E,z)e(z-z d 2 ) 



+M(E,z)Q( 



Z — Zn 



+ sm 2 e Q Ari(E,z)Q(zl-z) 
= cos 2 0< S) [MS t {E,z)+M! m {E,z) 
x {&(z - 4) + Q(z - z d )} 
+ S m 2 e Q ^ x (E,z)Q(z*-z), 

where G is the step function. On the other hand, if the relevant decay mode is Vi(pi 
+ X (helicity flips), the expected spectrum becomes 



(6) 



cos 2 6 Q Af Dl (E, z) + sin 2 8 G Afp 2 (E, z) 

\Ar° l (E,z)+Ar° 2 (E,z)e(z-z d 2 ) 



= COS 2 #0 



+N? 3 (E,z)Q(z-z a 3 
+ sin 2 9 Q Af? 2 



-z) 



= cos 2 9 { 



An 2 9 Q Af? 2 (E,z)G(z d 2 

Nl{E,z)+Nl{E,z) 
^-„- v --z d ) + {l-P H )Q{z-z d )} 
+ Af° x (E,z){(l-P H )&(z-z d ) 

+p H e(z-z d )}' 

+ sm 2 6 Q J^ x (E,z)G(z d -z). 
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Comparing Eqs. © and (J7J) with Eq. (J2J), the SRN flux with the neutrino decay is expected 
to be very different from the case of the pure neutrino oscillation. In particular, for the 



mode Vi 



+ X [Eq. ©], the SRN spectrum is expected to be hard since it contains a 



fair amount of u x . In the case of Ui — > Vj + X mode, the v e spectrum is also included, and 
then the corresponding upper limit is not as strong as that for z/j — > 9j + X mode. In the 
next section, we only consider the model, whose upper limit is the most severe among all 
the models considered, i.e., the decay mode z/j — > Vj + X. 

In general, the decay redshifts depend on the neutrino energy since the lifetime at the 
laboratory frame relates to that at the neutrino rest frame r via a simple relation 



Ti a b = Er/m. However, we believe that it does not make sense in discussing this point 
strictly, because the estimation of the SRN flux contains many other uncertainties as shown 
in Section IV Bl In order to obtain the lifetime of the mass eigenstates z/j from the decay 
redshifts, the typical neutrino energy E — 10 MeV is assumed with the following formulation: 



where the Hubble constant Hq is taken to be 70 km s™ 1 Mpc -1 , and A-dominated cosmology 
is assumed (Q m = 0.3, fl\ = 0.7). Since the exact value of the neutrino mass m is not 
known, the relevant quantities is the neutrino lifetime divided by its mass, r/m. 

IV. RESULTS 

Figure d shows the SRN flux for various parameter sets of decay redshifts (z$, zf) as a 
function of neutrino energy, which is calculated using Eqs. (0) and (jBJ). The solid curve in 
Fig. Q shows the SRN flux in the case that (zf, zf) = (5.0, 5.0), which represents the same 
flux as that without the neutrino decay; this SRN flux almost coincides with one obtained 
by the AST calculation. Then, we included the decay of the heaviest mass eigenstate u 3 by 
reducing the value of z%, with keeping u 2 stable. When z^ = 1.0 (dotted curve), the SRN 
flux at low energy region (E Pe < 35 MeV) deviates from the pure oscillation model (5.0, 5.0). 
This is because the neutrinos from supernovae at redshift larger than zf = 1.0 are affected 
by the z/3 — > v\ decay and it results in the increase of v e . Since the neutrino energies are 
redshifted by a factor of (1 + z)^ 1 owing to an expansion of the universe, the decay effect 
can be seen at low energy alone. When the value of zf is reduced to 10 -2 , the neutrinos 
even from the nearby sources are influenced by the ^3 — > vi decay, resulting in the deviation 
over the entire energy range as shown by the long-dashed curve in Fig. ^ If we add the 
z/ 2 — > v\ decay, it further enhances the SRN flux. 

In Table HJ we summarize the SRN flux integrated over the energy range of E Pe > 19.3 
MeV, for the each decay model. In the second column we placed the lifetime-to-mass ratio 
which corresponds to each decay redshift, which is obtained using Eq. (jHj). The corre- 
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(8) 



^(1 + n m z) (1 + zf - n A (2z + z 2 ) ' 
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10 20 30 40 50 60 

Neutrino Energy [MeV] 

FIG. 1: The SRN flux for various parameter sets of decay redshifts. Each label represents (zf, zf). 



TABLE I: The predicted SRN flux for various decay models and the corresponding SK limit (90% 
C.L.). Integrated energy range is E„ e > 19.3 MeV. The ratio of the prediction and the limit is 
shown in the fifth column. 



Model (z%,z$) (T 2 /m,T 3 /m) [s/eV] Predicted flux SK limit (90% C.L.) Prediction/Limit 



(5.0,5.0) 


(3.9 x 10 10 ,3.9 x 10 10 ) 0.43 cm" 2 s" 1 


< 1.2 cm" 


" 2 s- 


-1 


0.35 


(5.0,1.0) 


(3.9 x 10 10 ,2.4 x 10 10 ) 0.55 cm" 2 s" 1 


< 1.3 cm" 


" 2 s- 


-1 


0.42 


(5.0,0.2) 


(3.9 x 10 10 ,7.7 x 10 9 ) 0.93 cm" 2 s" 1 


< 1.2 cm" 


" 2 s- 


-1 


0.75 


(5.0, KT 2 ) 


(3.9 x 10 10 ,4.4 x 10 8 ) 1.0 cm" 2 s" 1 


< 1.2 cm" 


" 2 S" 


-1 


0.88 


(10~ 2 , 10- 2 ) 


(4.4 x 10 s , 4.4 x 10 s ) 1.4 cm' 2 s' 1 


< 1.2 cm" 


" 2 s- 


-1 


1.2 
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sponding 90% C.L. upper limit given by the SK observation at E Pe > 19.3 MeV [8| and 
prediction-to-limit ratio are also shown in the fourth and fifth columns, respectively. The 
predicted flux becomes larger as the decay is included, while the observational upper limit 
remains unchanged. At first sight, it is expected that the model (10 -2 , 10~ 2 ) is already 
excluded by the current observational data. 



n i i i r 



0.6 v... 



i r 





FIG. 2: A contour map of a prediction-to-limit ratio of the SRN flux which is projected against 
the lifetime-to-mass ratios (r^lvn^T^lm). 

Figure El shows a contour plot for a ratio of the predicted flux to the observational upper 
limit, which is projected against the lifetime-to-mass ratios (r2/m, r^jm). The area below the 
solid curve labeled as 1.0 is considered to be an excluded parameter region. We can confirm 
that our approach is very powerful to obtain the constraint on the neutrino lifetimes, because 
the best observed lower limit thus far (•> 10~ 4 s/eV) is much smaller than that shown in Fig. 
HI 10 10 s/eV), although there are still many uncertainties in this method as discussed in 
the next section. 
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V. DISCUSSION 



A. Other decaying modes 

Our discussions until this point were concerned with one specific decaying mode, z/ 3 — ► 

i>2 — > V\i where the daughter neutrinos carry nearly the full energy of their parent. We 
consider the other possible decaying models; first P 3 — > ^2,^2 —> &i is discussed. Since u 2 
state contains v e state by less fraction than U\, this decay mode gives smaller SRN flux, 
resulting in weaker upper limit in general. However, when the lifetime of one mode is much 
longer (shorter) than the other, i.e., z% <C zf or z 1 ^ ^> z%, the discussions for our reference 
decay modes i> 3 — > v 2 — »■ v\ are basically applicable to this case. 

When the daughter neutrino energy is considerably degraded, which is actually the case 
when the neutrino masses are strongly hierarchical, the observational upper limit for the 
each model is not as strong as the previous limit shown in Table |l] This is because the 
energetically degraded daughter neutrinos soften the SRN spectrum. 

Finally we consider the case that the daughter neutrinos are sterile species which does 
not interact with matter. If the lifetime of the mode is sufficiently short, then the obtained 
spectrum from each supernova can be expressed by 



which is smaller than the normal oscillation expression (J2J). Thus, also in this case, the 
observational upper limit will be looser. In consequence, all the other possible decaying 
models give weaker upper limit compared to our reference model. 

B. Supernova and supernova rate model uncertainties 

Again we restrict our discussion to our standard decay mode, and discuss whether the 
parameter region (r 2 /m, r 3 /m) ^ (10 10 , 10 10 ) [s/eV] is really ruled out by the current obser- 
vational data, as shown in Fig. |21 

The theoretical calculation of the SRN flux contains many uncertainties such as the 
supernova rate as a function of redshift z as well as the original neutrino spectrum emitted 
by each supernova. As for the supernova rate, since we have inferred it from the rest-frame 
UV observation, it may be severely affected by dust extinction, which is quite difficult to 




(9) 
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estimate. Thus for example, if we use the supernova rate model which is larger by a factor 
of 2 (even if this is actually the case, it is not surprising at all), then almost all the region 
of lifetime-to-mass ratio (V 2 /m, r 3 /m) will be excluded as shown in Fig. |21 

The uncertainty concerning the original neutrino spectrum also gives very large model 
dependence of the SRN flux calculation. Although we adopted the result of the Livermore 
simulation, as we have noted above it lacks the relevant neutrino processes, which reduce the 
average energy of v x to the value close to that of v e . If this is the case, the SRN spectrum as 
a result of the neutrino decay becomes softer compared to that obtained with the Livermore 
original spectra, which leads to weaker upper limit. In actual, we have calculated the SRN 
flux for various values of (z%, zf) assuming the original v x spectrum is the same as that 
of P e , for which we have used the Livermore data. As a result of the calculation for this 
extreme case, the obtained prediction-to-limit ratio is at most ~ 0.35, even when both of 
decay redshifts are sufficiently small. In consequence, considering the uncertainties which is 
included in the supernova rate or the neutrino spectrum models, we cannot conclude that 
the parameter region (T2/m,r 3 /m) ^ (10 10 , 10 10 ) [s/eV] is already excluded, although Fig. 
El indicates it. 



C. Future prospects 

Now, we consider the future possibility to set severer constraint on the neutrino decay 
models. The largest background against the SRN detection at SK is so-called invisible muon 
decay products. This event is illustrated as follows. The atmospheric neutrinos produce 
muons by interaction with the nucleons (both free and bound) in the fiducial volume. If these 
muons are produced with energies below Cherenkov radiation threshold (kinetic energy less 
than 53 MeV), then they will not be detected ("invisible muons"), but their decay-produced 
electrons will be. Since the muon decay signal will mimic the v e p —* e + n processes in SK, 
it is difficult to distinguish SRN from these events. Recent SK limits are obtained by the 
analysis including this invisible muon background. 

In the near future, however, it should be plausible to distinguish the invisible muon signals 
from the SRN signals; two different methods are currently proposed for the SK detector. One 
is to use the gamma rays emitted from nuclei which interacted with atmospheric neutrinos 
j^j]. If gamma ray events, whose energies are about 5-10 MeV, can be detected before 
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invisible muon events by muon lifetime, we can subtract them from the candidates of SRN 
signals. In that case, the upper limit would be much lower (by factor ~ 3) when the current 
data of 1,496 days are reanalyzed Another proposal is to detect signals of neutrons, 
which are produced through the v e p — > e + n reaction, in the SK detector. Actual candidate 

261 ] . Additional neutron signals 



for tagging neutrons is gadolinium solved into pure water 
can be used as delayed coincidence to considerably reduce background events. In addition, 
future projects such as the Hyper-Kamiokande detector is expected to greatly improve our 
knowledge of the SRN spectral shape as well as its flux. If a number of data were actually 
acquired, the most general and model-independent discussions concerning the neutrino decay 
would be accessible. Finally, we again stress that a great advantage to use the SRN is that 
the cosmologically long lifetime can be probed in principle. 



VI. CONCLUSIONS 



We obtained the current lower limit to the neutrino lifetime-to-mass ratio using the SRN 
observation at SK. Since the available data are only the upper limit on the flux integrated 
over E Pe > 19.3 MeV, the decay model on which we can set constraints is quite restricted, 
i.e., the decay models that gives the SRN flux which is comparable to or larger than the 
corresponding upper limit. 

Therefore, our reference decay model in this Letter is the two-body decay \>i — > i?j + X, 
with the daughter neutrinos which are active species and carry nearly full energy of their 
parent. The neutrino mass hierarchy is assumed to be normal (mi < m 3 ). The SRN 
calculation is also very sensitive to the adopted supernova model, i.e., neutrino spectrum 
from each supernova explosion; we adopted the result of the numerical simulations of the 
Lawrence Livermore group. Although their calculation is recently criticized since it lacks 
several relevant neutrino processes, we adopt it as a reference model. This is because it 
is only the successful model of the supernova explosion and we definitely need fully time- 
integrated neutrino spectra. 

Our calculations with these models shows that the neutrino decay model with 
(rz/m, Ts/m) ^5 (10 10 , 10 10 ) [s/eV] apparently gives the SRN flux that is larger than the 
current upper limit (see Fig. |2J). Since this value 10 10 s/eV is much larger than that of the 
limit obtained by the solar neutrino observation > 10~ 4 s/eV, our approach is shown to be 
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very powerful to obtain the implications for the neutrino decay. At present, however, owing 
to the large uncertainties such as supernova models, this lower limit 10 10 s/eV cannot be 
accepted without any doubt. Future experiments with the updated detectors and reanalysis 
with the refined method is expected to give us greatly improved information on the SRN 
flux and spectrum. In that case, the most general model-independent discussions concerning 
the neutrino decay would be possible. 
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